Nature of Phase Transitions of Superconducting Wire Networks in a Magnetic Field 
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We study I-V characteristics of periodic square Nb wire networks as a function of temperature 
in a transverse magnetic field, with a focus on three fillings 2/5, 1/2, and 0.618 that represent very 
different levels of incommensurability. For all three fillings, a scaling behavior of I-V characteristics 
is found, suggesting a finite temperature continuous superconducting phase transition. The low- 
temperature I-V characteristics are found to have an exponential form, indicative of the domain-wall 
excitations. 
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The presence of a quenched symmetry-breaking field 
is known to have important consequences on the order- 
ing of low temperature phases in many physical systems. 
A striking example is the pinning of a two-dimensional 
(2D) vortex lattice by a periodic potential in a super- 
conducting wire network. Without pinning, a 2D elastic 
vortex lattice in homogeneous superconducting thin films 
would not have long-range translational order at any fi- 
nite temperature [1] and cannot have long-range super- 
conducting phase coherence even at T = 0. A periodic 
pinning potential, however, when it is commensurate to 
the vortex lattice, can induce a gap in the low-energy ex- 
citation spectrum [2] and create a new thermodynamic 
phase [3] of a pinned 2D solid with true long-range trans- 
lational order, and with superconducting phase coherence 
[4] . In the presence of a high-order commensurate (or in- 
commensurate) potential, the competition between the 
vortex-vortex interactions and the vortex-network inter- 
actions leads to a whole new class of problems [5-7] . For 
example, the vortex lattice may: (a) become a 2D "fioat- 
ing" soHd and again lose its translational order and su- 
perconducting phase coherence at any temperature [5]; 
(b) form a metastablo "glassy" phase [6]; (c) be pinned 
in a commensurate phase [7] and thus superconducting 
at low temperatures. The issue is far from being settled. 

Closely related to the low-temperature thermodynamic 
phase, the nature of the superconducting transition of a 
superconducting network in a magnetic field is also not 
well understood. At filling / = 1/2, where / is the frac- 
tion of a flux quantum 4>q = hc/2e per plaquette, for 
example, the vortex conflguration of the ground state of 
the system has a checker-board pattern [4]. Thus, the 
ground state of the system has the discrete symmetry of 
the two-fold degeneracy as well as the continuous sym- 
metry of an arbitrary global phase change. Two types 
of excitations are possible: vortex-antivortex pairs of the 
continuous phase variable and domain walls of the two 
ground states. If the two types of excitations do not inter- 
act, they should lead to two independent transitions [9]: 
a Kosterlitz-Thouless (KT) transition [8] in the underly- 



ing network and an Ising melting transition in the vortex 
lattice. The superconducting transition will be deter- 
mined by the one with the lower transition-temperature. 
Interesting physics arises when the two types of excita- 
tions do couple [10]. Due to the screening of the vortex 
interactions by the domain walls, when the domain wall 
energy goes to zero at an Ising-like transition, a sponta- 
neous generation of domain walls will unbind the vortex 
pairs and induce a KT-like transition [10]. In this sce- 
nario there will be a single superconducting transition 
of a new universality class [10,11]. In spite of intense 
theoretical effort on this problem, the exact nature of 
this transition remains unknown [11]. We show in this 
paper that the superconducting transition at f = 1/2 is 
continuous and is strongly influenced by the Ising-like ex- 
citations (domain walls), but is not an Ising transition. 
Moreover, we flnd that the phase transitions at higher 
order commensurate or incommensurate fillings such as 
/ = 2/5 and f = 0.618 are, surprisingly, also similar in 
nature. 

We describe here an experiment on a periodic square 
superconducting wire network. The network sample in 
our experiment is fabricated from high-quality sputtered 
Nb film. The microfabrication consists of standard steps, 
i.e., electron beam lithography, lift-off and reactive ion 
etch. The final sample has 800 x 800 cells, with a lat- 
tice constant 1 jum, wire width 0.25 //m and thickness 
0.1 /im. SEM (scanning electron microscope) examina- 
tions of the final sample show that the overall periodicity 
is maintained very well, as shown in Fig. 1. The zero- 
field transition has (at 0.5R„)k 8.920 K and width 
((10% - 90%)i?„) « 10 mK. The normal-state resistance 
Rn « 1.2 f2/square. The standard four-probe technique 
is used for the transport measurements. To ensure cur- 
rent uniformity, we use large contact pads (silver paint), 
each covering the entire edge of the network. The volt- 
age contacts are lithographically patterned with the net- 
work, both are 5 fim in width and are of the same ma- 
terial Nb. The close-loop temperature regulation is done 
by a Linear Research LR-130 controller and the voltage 
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is measured by a transformer-coupled PAR- 124 A lock-in 
amplifier with a square-wave current at 10.2 Hz. To test 
possible heating effects on the I-V curves, we use differ- 
ent duty cycles and find no measurable differences. The 
temperature reading accuracy around 9 K is about 1 mK, 
although the temperature stability is better than 0.5 mK 
(determined by using the sample itself as a thermome- 
ter). 




FIG. 1. SEM micrographs showing the structure of the su- 
perconducting Nb wire network used in this work. 

Fig. 2 shows the sample resistance as a function of mag- 
netic field at various fixed temperatures. As the sample 
is cooled from the normal state, the sample resistance 
first dips at / = and other integer fillings: upon fur- 
ther cooling, the resistance dips at / = 1/2 and other 
fillings. Taking advantage of the variation of the sample 
resistance with /, one can determine a mean-field phase 
boundary [12] that resembles some features of the inter- 
esting Hofstadter spectrum [13]. That phase boundary 
is not, however, the true phase boundary of the super- 
conducting transitions due to the neglect of the vortex 
fluctuations. 
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FIG. 2. Sample resistance R{= V/I) as a function of 
filling factor / at various temperatures: from top, 8.920, 
8.912, 8.907, 8.901, 8.895 and 8.890 K. The test current was 
lac = 10/iA at 10.2 Hz. 

Since the emergence of long-range superconducting 
phase coherence is characterized by nonhnear I-V charac- 
teristics, we study how the I-V curves behave as the sam- 
ple is cooled in a flxed fleld. As a calibration, we measure 
the temperature dependence of the I-V characteristics at 
/ = 0, 1, and 3 [14]. We find that the I-V curve changes 
from ohmic (at low current) to power-law with decreas- 
ing temperature. The low-temperature power-law I-V 
curves are characteristic of the vortex-antivortex pairs 
bound by a logarithmic interaction. The temperature 
dependence of the power-law exponent a {V ^ I"^) shows 
the familiar jump at a « 3 [15]. Other aspects of the 
data for / = 0, 1, and 3 are also consistent with earlier 
works [15]. We thus conclude that at / = 0, 1, and 3 our 
system undergoes a KT transition. 

In sharp contrast to the power-law I-V curves at / = 0, 
1, and 3, the low-temperature I-V characteristics at non- 
integer fillings are found to have an exponential form and 
the I-V curves behave similarly for different non-integer 
fillings as a function of temperature. Fig. 3 shows the 
I-V curves for f = 2/5, 1/2, 0.618, and 5/2 [14]. At high 
temperatures and for all these filhngs, the I-V curves 
are ohmic at low current and concave upward at high 
current. At low temperatures, the I-V curves become 
progressively more concave downward. The high-T and 
low-T behaviors are separated by a F ~ P power-law. 
The lowest temperature curves for all four fillings in Fig. 
3 can be fitted with a form V ~ /oxp(— Ix//), where It 
is a fitting constant. This exponential form can be de- 
scribed by a simple domain- wah nucleation process [16]: 
assuming the ground state vortex lattice is pinned by 
the underlying network, the motion of the lattice is via 
a thermally activated domain nucleation process. The 
energy cost for creating a domain scales as crL, where a 
is the wall energy per unit length and L the length of 
the domain wall; the energy gain from transport current 
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/ scales as IL^. Thus the energy barrier for growing an 
unbounded domain scales as and this process gives 
an exponential form V ~ Iexp{—lT/I) for I-V curves. 
As discussed earlier, due to the symmetry of the ground 
states for non- integer fillings, both vortex pairs and do- 
main walls are possible excitations. Note that the energy 
cost to creat free vortices from bound pairs diverges loga- 
rithmically with vanishing current, while the cost to grow 
an unbounded domain diverges algebraically. Therefore, 
if the two types of excitations are independent the low- 
temperature I-V curves should be power-laws. The ob- 
served exponential I-V curves are then cither due to the 
coupling between the two types of excitations, or a result 
of a vanishing wall energy cr — > close to the transition. 
For / = 1/2, we note that an early experiment [17] on 
Josephson arrays reported power-law I-V curves at low 
temperatures, and a recent experiment [18] on Al wire 
networks showed I-V curves that appear to be exponen- 
tial. The origin of this difference between the two systems 
is not clear. 
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FIG. 3. Temperature dependent current-voltage character- 
istics at various fillings / =2/5, 1/2, 0.618, and 5/2. The I-V 
curves correspond to the following temperature ranges: for 
/ = 2/5, T =8.900 K (top I-V) to 8.850 K (bottom I-V); for 
/ = 1/2, T =8.880 K to 8.847 K; for / = 0.618, T = 8.890 K 
to 8.845 K; for / = 5/2, T = 8.905 K to 8.855 K. The dashed 
lines are hand-drawn to indicate the V ~ power-law. 

The exponential form of I-V characteristics suggests 
superconducting phase coherence at low temperatures in 
the wire network. For f = 1/2 and 2/5, this may be un- 
derstood as the long-range ordering in the pinned ground 
state vortex lattice. For higher-order commensurate or 
incommensurate fillings, e.g. / = 0.618, the existence of 
an ordered ground state could be an issue [6] . For adsor- 
bate systems [7] a strong pinning potential can force an 
incommensurate lattice into registry with the substrate 
potential. The ground state is a pinned commensurate 
lattice [7]. However, in our system the vortex density is 
fixed by an external field and the screening is very weak 
in the temperature range of the experiment. The ground 



state of the system may be cither commensurate domains 
separated by domain walls, or ordered quasiperiodically. 
The exponential I-V curves found in this work and in a 
previous one [19] suggest a true superconducting state at 
low temperature for / = 0.618. This would imply that 
the vortex lattice is pinned at low temperature, either 
by the network or by the network irregularities which 
undoubtedly exist in our sample (a worst-case estimate 
from the SEM photo in Fig. 1 gives ~ 10% in areal fiuc- 
tuations among the cells ) and that of the previous work 
[19]. 



The I-V data in Fig. 3 suggest phase transitions from 
the high-temperature resistive state to low-temperature 
superconducting state. If the superconducting transi- 
tions are continuous, one would expect to sec scaling 
behavior in I-V characteristics [20,21]. In the follow- 
ing we attempt a scaling analysis of the I-V data in 
Fig. 3. Approaching the transition, the phase correla- 
tion length ^ (not to be confused with the bulk super- 
conducting coherence length) in the network diverges as 
^ ~ |T — Tel"". Dimensionality arguments [21] suggest 
a scaling relation for the I-V (or j-E) curves, in 2D, as 
{E/j)\T - T,]-^" = G±{j\T - T,|-''), where z is the dy- 
namical critical exponent and G± the scaling functions 
above and below T^. At T = T, and in 2D, V ~ 
thus we can determine z directly from the I-V curves 
in Fig. 3. For all the cases which we studied carefully, 
z = 2; this includes f = 2/5, 1/2, 0.618, and 5/2. In 
the scaling analyses, we fix the parameter z = 2, as de- 
termined independently from the power law, and keep ly 
adjustable until a good collapse of data is achieved. The 
scaled data are shown in Fig. 4. The four sets of I-V 
curves in Fig. 3 are found to show good scahng behavior, 
which strongly suggests a finite-temperature continuous 
superconducting transition for all these fillings. This is 
surprising because one would expect very different types 
of transitions for these different fillings. 
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FIG. 4. Scaling plots for / =1/2, 2/5, 0.618, and 5/2. The 
insets showing the values of z, v and Tc used to scale the data 
(see text). 
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At / = 1/2, one may have a continuous transition 
[9-11]. But for / = 2/5, one expects [22] a strongly 
first-order melting transition. At / = 0.618, a glass-like 
freezing into a metastable state may occur at a finite 
temperature, but no finite-temperature phase transition, 
thus no scaling behavior, is expected [6], in contrast to 
the interpretation in a previous work [19]. As discussed 
earlier, a true superconducting phase at low tempera- 
tures for / = 0.618 may be a result of the network ir- 
regularities. The fact that there is no pronounced dip 
feature at / = 2/5 on the R vs. / curves (Fig. 2) sug- 
gests that the / = 2/5 ground state may also be affected 
by such disorder. Since the first-order transitions can 
be driven continuous by quenched disorder [23,24], espe- 
cially in 2D [23], it is possible that the continuous tran- 
sitions for / = 2/5 [25] and 0.618 are due to the areal 
disorder in the network. For / = 1/2, we note that the 
critical exponent « 1.7 is different from the 2D Ising 
value {v = 1) and from the nonuniversal value {v « 0.8) 
of a coupled XF- Ising model [11]. In random-bond 2D 
8-state Potts model, the transition was found [24] to be 
driven by disorder to the Ising universality class. If weak 
disorder plays a role in driving a first-order transition 
continuous in our network, as conjectured above, an in- 
teresting question is whether the transitions at / = 2/5, 
1/2, and 0.618 are in the same uniA-ersality class. Unfor- 
tunately, due to a systematic error (~ 4 mK) in locating 
Tc, the error bar for exponent f can be 0.3 to 0.5. Thus 
it is not clear whether the variation in i/ for different fill- 
ings represents a real change in their universality classes, 
or they are in fact in a single universality class of cither 
/ = 1/2 or one controlled by disorder. In any case, the 
issue of disorder-induced continuous transition [23,24] in 
a 2D network deserves more careful studies and is cur- 
rently under investigation. 

Another remarkable feature of the data is the F ~ 7^ 
power-law for / = 2/5, 1/2, and 0.618 at the transition, 
which is typical of a KT transition, but which occurs 
while the I-V curves clearly suggest domain-wall exci- 
tations just below the transition. As mentioned earher, 
the domain walls have the effect of screening the vortex- 
antivortex interactions, it is possible that an Ising-like 
transition at which the domain wall energy goes to zero 
may trigger a KT transition [10]. It is not clear, however, 
why the universal power of 3 which is related [26] to the 
universal jump in the helicity modulus of a KT transition 
should appear at such a transition. 
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